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LONG-TERM MOTION OF A LUNAR SATELLITE 


I. INTRODUCTION 

The purpose of this report is to examine two special methods of studying 
the long-term motion of a lunar orbiter. The Hamiltonian used is the long- 
period Hamiltonian (in which the argument of pericenter and inclination of the 
satellite appear but other angular variables have been eliminated) as derived 
by Kozai and Giacaglia, et. al., (References 2, 3). It includes the perturbation 
of the earth as a point mass and the principal part of the oblateness of the moon. 

This long-period Hamiltonian is 




l 


F = — KjL (5 - 3^2) 


= -C, 



K 2 v 3 




+— K.L(1 
2 




2g 


(1) 


where 

r/ 2 = 1 - e 2 

e = eccentricity of the satellite's orbit 
L = //u. a 
G = L /l - e 2 
H = G cos i 

/j. = gravitational constant times the mass of the moon 
= 3.6601891 x 10“ 3 decamegameters 3 /centiday 2 

a = semimajor axis of the satellite's orbit (decamegameters) 

i = inclination of satellite orbit plane to the moon' s equatorial plane 

g = argument of pericenter of the satellite 
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Kj = (centiday -1 ) 

= mean motion of the moon 
= 2.2802713 x 10" 3 (radians/centiday) 

n = m 2 L -3 = mean motion of the satellite 

e = ratio of sum of masses of earth and moon to the mass of the earth 
= 1.0123001 

K, = — J,b 2 n 2 (decamegameters 2 /centiday 2 ) 
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J 2 = principal part of the oblateness of the moon 
= 2.41 x 10 " 4 (Reference 2) 

b = mean radius of the moon 
= 0.1738 decamegameters 


Since the Hamiltonian (1) is not an explicit function of time, it has a constant 
value, F = 1/6 C . 

Special methods are needed for the problem of the long-term motion since 
the standard procedures — von Zeipel's method and the method of successive 
approximations — break down at this point for the lunar satellite (e.g. see Ref- 
erence 2). The special methods, which involve harmonic analysis and elliptic 
integrals, to be considered in this report were suggested many years ago by 
E.W. Brown in his paper "On the Stellar Problem of Three Bodies" (Reference 1). 

In the method involving harmonic analysis it is assumed that 17 2 and dt/dg 
may be represented by cosine series with the argument of pericenter g as the 
variable. Integrating dt/dg gives g as a sum of a linear expression of time 
and a sine series in g . On reversing this expression we have g as a function 
of time. Then tj 2 is approximated as a cosine series in time. In the second of 
these methods proposed by Brown, 


. 1 3F . BF 

77 = , g = , 

L Bg BG 


2 
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and an approximation of cos 2g found by neglecting the oblateness of the moon 
are used to give d77 2 /dt and dg/drj 2 as functions of v 2 alone. These result in 
elliptic integrals of the first and third kinds as the solutions for rj 2 as a function 
of time and g as a function of rj 2 (and therefore of time) . 

For the harmonic analysis method it is necessary to solve the Hamiltonian 
for v 2 corresponding to chosen values of the argument of pericenter. To be 
able to compare the values of rj 2 and g as obtained from harmonic analysis and 
elliptic integrals the portion of the Hamiltonian containing the oblateness of the 
moon was neglected, both in solving for rf and evaluating dt/'dg. Setting J 2 = 0 
results in the following quadratic equation in x = rj 2 : 


3(1-5 co s 2g) x 2 - [5 + 9v 2 -15cos2g(l + v 2 ) + C'] x + 1 5 v 2 (1 - cos 2g) = 0 


where 


C' = 2(5 - 6 r} 2 + 3 v 2 ) 


and v = H/L are constants, rj 2 is the value of rj 2 when g = 0. If the oblate- 
ness of the moon is included in the Hamiltonian, however, an approximation used 
in this report is to represent (1 - e 2 )~ 3/2 by its Maclaurin series terminating 
the series to give the quadratic equation in e 2 (for small to moderate values 
of e ) 


ae 4 + fie 2 + y - 0, 


where 


a = 3KjL + — K (-1 +15v 2 ) - 15K 1 Lcos2g 
4 

fi = (KjL + K 2 ) (-1 + 9^ 2 ) + C + 15 KjL (1 - v 2 ) cos 2g 
y- 2(KjL + K 2 ) (-1 + 3v 2 ) -C. 
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This quadratic is then solved for r, 2 - \ - e 2 , given specific values of cos 2g 


V 


2 




(the sign is chosen so that rj 2 lies between zero and one) . 


II. THE HARMONIC ANALYSIS METHOD 


Given the initial conditions 

L = 0.06 
G = 0.055 
H = 0.05 
g = 0 

we find 

e = 0.39965264 

a = 0.98 355573 decamegameters 
i = 24.619974° 

Kj = 3.105573 x 10 ' 5 
K 2 = 2.1003149 x 10- 8 
C = 7.6893300 x 10' 6 - 


We assume 


V 


2 


L 

k > 0 


a^, cos 2kg 


and 


dt 

dg 


L 

k >0 


terminating the series with co s 8 g so that 


V 2 - a 0 + a x cos 2g + a 2 cos 4g + a 3 cos 6g + a 4 

— = b 0 + b x cos 2g + b 2 cos 4g + b 3 cos 6g + b 
dg 


b k cos 2kg, 

cos 8g 
cos 8 g. 


4 


4 


Solving the Hamiltonian for rj 1 given the five values 2g =0, 77,77/2, 77/3,277/3 , 
we have the following results: 


2g 

0 

77 

77 

2 

77 

3 

277 

3 

V 2 , K 2 = 0 

.84027778 

.77996317 

.80665828 

.82299779 

.79219924 

V 2 , k 2 4 0 

.84027778 

.78055119 

.80702653 

.82317074 

.79270692 


Corresponding to each of the five special values of 2g with the associated values 
of v 2 , dt/dg is evaluated using 


dg 

dt 


3F 

Kj77 

BG 


cos 2g 



Then the following table may be used to solve for the coefficients in the series 
for dt/dg : 


2g 

0 

77 

77 

2 

77 

"3 

2 77 
3 


8781.8656 

6732.7808 

8268.1582 

8735.5882 

7538.4678 


8663.7170 

6662.5543 

8155.7627 

8611.8994 

7446.8692 


Letting f k denote the k th value of rj 1 in the table, the equations 


f 

f 


1 

2 


a 0+ 3 1 + a 2 + a 3 


a 0 3 1 + a 2 a 3 


+ 3 

+ a 


f 3 = *0 
f 4 = a 0 


- a 


2 


+ a 


4 



- a. 


+ a. 
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1 
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may be solved for the coefficients a k in the terminated series expansion for v 2 . 
For K 2 = 0 we then have 

a„ =|< f l + f 2 >+ + f s> = 80843917 
a, =!< f . - <2 * f 4 - f 5 > = .03637106 

a 2 =-j(f, +f 2 ) -| , 4 = 00173110 

a = i(f,-f 2 )-I(f 4 -f 5 ) = - 00021374 
6 6 ^ 

a 4 = f 3 - a 0 + a 2 = - .00004979. 

Thus, if K 2 = 0, 

rj 2 - .80843917 + .03037106 cos2g + .00173110 cos 4g (2) 

- .00021374 cos 6g - .00004979 cos 8g 

while, if K 2 ^ 0, 

T] 2 - .80876405 + .03006347 cos 2g + .00169 398 cos 4g 

- .00020018 cos 6g - .00004354 cos 8g. (2') 

Similarly, the five coefficients b fc in the terminated cosine series for dt/dg may 
be evaluated, so that if K 2 = 0, 

ll = 8010.4597 + 1082.0684 cos 2g - 255.4175 cos 4g 
dg 

- 57.5260 cos 6g + 2.2810 cos 8g 
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and, if K 2 ^ 0, 


H 


_ = 7907-3015 + 1055.3976 cos 2g - 246.3136 cos 4g 
dg 

- 54.8163 cos 6g + 2.1476 cos 8g. 


Integrating this expression for dt/dg gives 

b l b 2 b 3 b 4 

t = b n g + sin 2g + — sin 4g + sin 6g + — sin 8g 

0 2 4 6 8 

(the constant of integration is zero in this example because of the initial condi- 
tion g = 0) . Then 


1 b i b 2 b 3 b 4 

g- — t - sin 2g - sin 4g sin 6g sin 8g. 

b„ 2b„ 4b„ 6b„ 8b 0 


In particular, 


(K 2 = 0) g = .12483678 x 10' 3 t - .06754097 sin 2g 


+ .00797137 sin 4g + .00119689 sin 6g - .00003559 sin 8g 


(K 2 / 0) g = .12646539 x 10" 3 t - .06673563 sin 2g + .00778754 sin 4g 
+ .00115539 sin 6g - .00003395 sin 8g . 


Since we may write g = r + cp( g) , where r is a linear function of time and 
9(g) is a trigonometric function in g , Lagrange's expansion theorem may be 
applied to reverse the equation and give g as a function in time. By this 
theorem, if 


y = x + 9(y) , 


7 




then a series may be used to express y as a function of x: 


y = x + cp(x) + — — (cp 2 (x))+— (? 3 (x)> -t • • 

2 ! 3x 3 ! 


In this case, only the first four terms of this series were taken and all terms 
after sin 8g were eliminated. 

Hence, 


where 


g= T+ ( a+ 4 + t) sin 2t + ( b + I) sin 4 


+ ct I t - sin 6 t + fd + sin6 


k 2 = ° 


k 2 f 0 


-.067540965 

.0079713748 

.0011968934 

-.000035595648 


.06673563 

.00778754 

.00115539 

.00003395 


A = - 2(ab + be + cd) 
V = 2(a 2 - 2ac - 2bd) 


J =6a(b -d) 


D = 4(2ac + b 2 ) 


E - 6 J - — - b 2 a - c 2 a - d 2 a bed + + bad 


a 2 d b 3 c 2 d 

~2 2 2 ~ 
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4k 


M = 54 - b 2 c - bed — d 2 c - adb + 

N = 48 ^ba 2 - 2 abc - 2b 2 d - be 2 - 2 a 2 d - d 3 - 2dc 2 

Thus, for K. = 0, 

g = r - .06685092 sin 2r + .01255112 sin 4r 

- .00091832 sin 6 r + .00007192 sin 8 r ^ 3 ^ 

where r = 0.12483678 x 10' 3 t ( t in centidays), while if K 2 4 0, 

g - r - .06606959 sin 2 r + .01225829 sin 4 r - .00088524 sin 6t 

+ .00006836 sin 8r. (3’) 

( T = .12646539 x 10“ 3 t, t in centidays). 

Letting 2r take on the five special values 0, n, -tt/2 , n/3, 2v/3, these 
equations for g give corresponding values of 2g; from these we may evaluate r/ 2 
using the terminated cosine series (2) and (2’). This gives the following results: 


2 r 

0 

17 

77 

Y 

77 

Y 

277 

3 

g , K 2 = o 

0 

77 

*2 

.71946556 

.47651148 

.97849565 

V 2 , K 2 = 0 

.84027778 

.77996317 

.81080056 

.82570548 

.79556239 

g , K 2 4 0 

0 

77 

Y 

.72021381 

.47693763 

.97942282 

v 2 , K 2 4 o 

.84027778 

.78055119 

.81107413 

.82582292 

.79599495 
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Assuming r / 2 may be written as a cosine series in r and taking the first 
five terms of the series, we may evaluate the coefficients of these terms as 
above. Thus, we have for K 2 = 0, 

V 2 = .81046278 + .03015257 cos 2 r - .00034004 cos 4r 
+ .00000474 cos 6 r - 00000226 cos 8 r 

(r = .12483678 x 10' 3 t) (4) 


and for K 2 ^ 0, 


rf = .81074412 + .02985152 cos 2 r - .00032982 cos 4r 

+ .00001177 cos 6 r + .00000019 cos 8r. (4') 

(r = .12646539 x 10" 3 t) 


In this way both g and rj 1 may be found as functions of time. 

Tables 1 and 2 give the values of rf and g. Equations (4) and (4') are 
plotted against time in Figure 1. This shows the shift that occurs when the 
oblateness of the moon is considered in the long-period Hamiltonian. Figure 2 
shows similar results for g as given by Equations (3) and (3'). 

As may be seen from these figures, the effect of the moon's oblateness on a 
satellite at this distance from the moon (5.659 moon radii) is not great. For a 
satellite near the surface of the moon, the effect of the moon's oblateness on the 
eccentricity and argument of pericenter of the satellite would be much more 
pronounced. 
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Table 1 

V 2 , g Harmonic Analysis (K 2 = 0) 


t 

(days) 

V 2 

g 

Radians 

g 

Degrees 

F 

x 10" S 

F - C 
x 10 ' 10 

0 

.84027782 

0 

0 

.76086524 

0 

20 

.83675070 

.2273639 

13.02699 

.76086623 

.089812 

40 

.82692477 

.4543801 

26.03406 

.76086372 

-.160298 

60 

.81298783 

.6850423 

39.25003 

.76086474 

-.059685 

80 

.79215000 

.9287134 

53.21135 

.76086799 

.266027 

100 

.78627455 

1.1952296 

68.48161 

.76085606 

-.927684 

120 

.78029692 

1.4843190 

85.04521 

.76086390 

-.143245 

140 

.78191710 

1.7799106 

101.98136 

.76085884 

-.649720 

160 

.79068695 

2.0606739 

118.06792 

.76086051 

-.481463 

180 

.80424327 

2.3168499 

134.33726 

.76086867 

.334239 

200 

.81911010 

2.5537270 

146.31778 

.76086276 

-.2575 

220 

.83167979 

2.7820046 

159.39712 

.76086533 

0 

240 

.83906558 

3.0090020 

172.40312 

.76086577 

.043769 

260 

.83965478 

3.2365252 

185.43923 

.76086555 

.022168 

280 

.83332190 

3.4636120 

198.45035 

.76086267 

.051159 

300 

.82144059 

3.6913262 

211.49741 

.76086784 

-.266027 

320 

.80671501 

3.9260582 

224.94656 

.76086301 

.250679 

340 

.79268607 

4.1783800 

239.40354 

.76086301 

-.232489 

360 

.78291847 

4.4552836 

255.26894 

.76085703 

-.829345 

380 

.78002546 

4.7497567 

272.14101 

.76086506 

-.027284 

400 

.78480573 

5.0413928 

288.85053 

.76085565 

-.967474 

420 

.79594596 

5.3120934 

304.36053 

.76086641 

.107434 

440 

.81050099 

5.5589943 

318.50691 

.76086596 

.062527 

460 

.82482343 

5.7910522 

331.80284 

.76086308 

-.225099 

480 

.83551665 

6.0182398 

344.81974 

.76086616 

.082991 

500 

.84017939 

6.2454964 

357.84058 

.76086530 

-.003410 

520 

.83780875 

.1897622 

10.87257 

.76086611 

.077875 

540 

.82891087 

.4167136 

23.87593 

.76086439 

-.093791 
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Table 2 

ri 2 , g: Harmonic Analysis (K 2 ^ 0) 

t 2 g g F AF 

(days) v Radians Degrees x 10“ 5 x 10' 8 

0 .84027778 0 0 .76393302 0 

20 .83668256 .23048543 13.205842 .76898250 .0495106 

40 .82670981 .46072723 26.397726 .76906286 .1298758 

60 .81265809 .69497792 39.819302 .76911933 .1863440 

80 .79790366 .94291820 54.025233 .76917259 .2395950 

100 .78620654 1.2142525 69.571544 .76921679 .2837964 

120 .78072811 1.5074170 86.383185 .76926319 .3201989 

140 .78300817 1.8058868 103.46969 .76923517 .3021796 

160 .79240119 2.0877218 119.61765 .76919183 .2588421 

180 .80631778 2.3446295 134.33738 .76914220 .2092065 

200 .82111800 2.5832362 148.00853 .76908830 .1553075 

220 .83316153 2.8143461 161.25016 .76901856 .0855607 

240 .83963119 3.0445164 174.43794 .76894303 .0100385 

260 .83905608 3.2751774 187.65384 .76895149 .0184968 

280 .83156597 3.5053372 200.84103 .76903174 .0987483 

300 .81886887 3.7369943 214.11400 .76909695 .1639591 

320 .80395451 3.9777790 227.90995 .76915085 .2178580 

340 .79050779 4.2384158 242.84334 .76919868 .2656861 

360 .78212960 4.5239208 259.20157 .76924155 .3085574 

380 .78109380 4.8229269 276.33336 .76924995 .3169589 

400 .78771252 5.1136422 292.99012 .76920982 .2768274 

420 .80013334 5.3809400 308.30515 .76916465 .2316596 

440 .81502776 5.625948 322.34308 .76911089 .1778971 

460 .82865910 5.8590324 335.69783 .76905185 .1188595 

480 .83771752 6.0893168 348.89215 .76896972 .0367322 

500 .84018579 .03656091 2.094786 .76893447 .0014779 

520 .83548432 .26695049 15.295136 .76899601 .0630166 

540 .82467794 .49738946 28.498317 .76907304 .1400508 


2 











KJs z#ro, 






is zero. S; K^is not zero, 
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m. THE SOLUTION WITH ELLIPTIC INTEGRALS 

The second method suggested by Brown, starting with 77 and g, results in 
elliptic integral expressions for rj 2 and g, but for this solution the oblateness 
portion of the Hamiltonian is neglected. 

From the long-period Hamiltonian ( 1 ) and the definition of the Delaunay 
variables, we have 


£77 \ dG 1 BF_ n 

dt 'Ldt _ L3g 



sin 2g. 


Setting K 2 = 0 in the Hamiltonian gives 

W= (5- 3 rj 2 ) ^1 +±lj+ 15(1 -tj 2 ) ^1 -Z!jcos 2g, 


( 5 ) 


( 6 ) 


where the constant W is evaluated for g = 0. Then, w = 2(5 - 6r) 2 + 3v 2 ) where 
77 x 2 is the value of rj 2 when g = 0. On solving the Hamiltonian ( 6 ) for cos 2g, 
we have 


cos 2 g = 


15(1 - 77 2 ) 



(5 - 3t7 2 ) 



+ 2(5 - 6 V 2 +3v 2 ) 


which may be substituted in ( 6 ) forsin 2 g = [1 - cos 2 2 g] 1/2 to give 


-77 = 27777 = + 4^6 Kj [(x - Xj)(x - x 2 )(x - x 3 )] 1/2 . 


where x = rj 2 , x x = 77 ^ and the roots x 2 ,x 3 of the polynomial are defined by the 
equations 


5 2 

* 2*3 


Xj + x 3 =1 (5 + 5v 2 — 2 xj). 


In this expression for dx/dt the minus sign is used when sin 2g is positive 
and the plus sign when sin 2 g is negative. 
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If the roots of the polynomial are relabeled to give a strict ordering 
x | < x 2 < with the relationship xj < x < x' 2 , we have 

- dx — = + 4Kj /6 dt. 

/(X - x') (x - x' ) (x - x') 


Integrating both sides gives 


f 


dx 


A* - x i) ( x - x 2) ( x - x ' 3 ) 


= T 4K. /6 


£ 


dt 


where x Q is the initial value of x and t Q may be taken as zero. 
Set 


x = 


X. - X, y 


1 - y 2 


dx = 


2(x' - x ')y dy 


(i -y 2 ) 


2x2 


k = 


X 2 - X l 


x 3 - X t 


and we have 


T 4Kj t / 6 



dy 

/(I - k 2 y 2 ) (1 - y 2 ) 
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m 


or 





Let us define 


v = 2K, t 




Then 



17 


so that 


u + V 



dy 

Al - y 2 )(l - k 2 y 2 ) 


Thus we may solve for x as a function of time, using the Jacobi elliptic function 
of u and v: 



sn (u i v) 


x = Xj' + (x 2 ' - x x ') sn 2 (u ± v) 


= x/ + (x 2 ' - x/) 


sn u cn v dn v ± sn v cn u dn u" 
1 - k 2 sn 2 u sn 2 v 


For the given initial conditions we have 

x j = .77996317 
x' = .84027778 

x' = 1.48392569 
k 2 = .08567872 
sn u = 1 
cn u = 0 


x — .77996317 + .06031461 


cn v dn v 

1 - .08567872 sn 2 v 
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To evaluate g , the Hamiltonian is solved for cos 2g, giving 


cos 2g = 


15(1 - x) 1 


-t) 


-1 


c_ + i!i + <5 - 3x> 


_K,L K,L 

This expression is substituted for cos 2g in 


3/2 


(-¥)' 


g = * 


__L 

/ x 


(* -¥) - s(x - £) COS 2g] - A (l - 


<Bince 


4 = ii • = + 4/? K /QOO 

dx dt dx 


where Q(x) = (x - Xj )(x - x 2 ) (x - x 3 ) , this substitution gives 


/xQ(x) ' 

+ (5 - 3x) 


V ^ (x-i^ 

r c + 2r 2 

1 - 3v 2 x _1 

' AQ(x) \ x ' 

Lk x L k,l 

x 3/2 

.^iLd-x) fi 

..at. 

k 2 

x /JL \ 

* /j 

Lx 2 /Q(x) 


(•-¥) 


Integrating this expression, we have 


+ 4/6 K, g = — 


3 LL 


— + 2(1 - 3v 2 )K 


■]£ 


dx 


x 0 (1 - x) v xQ(x) 


_^[c 

" 3U 


+ 2(3v 2 - 5) K 


If <! 

l JJ. (X - C J ) 


dx 

/xQ( x) 
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_i 

3 


+ 2(1 + 3v 2 ) Kj 


■]f 


dx 

i/xQ(x) 


2 K 2 


r 

•'X- 


dx 


+ - — ( l-3v 2 ) 

3 L J x 0 (1 - X) /QOO 


_ ±*2 f 

3v 2 L J 


? K 9 

- -r o 2 - 3 y4 - 2) 

3 v 2 L 


s. 


dx 


^ x/Q(x) 


(x - v 1 ) v/Q( x) 


1 K 2 

— — (7 + 6 y2 ) 
3 L 


I 


dx , 2 K 
+ 3 yi 


x 0 X 2 /Q(x) 


H 


dx 


x J v Q(x) 


In actually evaluating g only the first three summands of the preceding were 
used since the oblateness of the moon was omitted from the Hamiltonian. The 
elliptic functions were evaluated using Reference 6. 

IV. CONCLUSIONS 

The results of these two special methods involving harmonic analysis and 
elliptic integrals are compared using the given initial conditions and the harmonic 
analysis solutions (3) and (4) for rj 2 and gwhen K 2 = 0. These results are tabulated 
in table 3 together with the differences between the values obtained by the two 
methods. 

In table 2 are the values of 17 2 and g when the oblateness of the moon is 
retained in the Hamiltonian and a Maclaurin series is used to give a quadratic 
equation in e 2 . The value of the Hamiltonian (1) was computed for 77 2 and K 
(at 20 day intervals) and compared to the initial value ( rj 2 = .84027778, g = 0) 

C = .76893300 x 10 ~ 5 . 


Although g is periodic in this example and all values are admissible, there 
is the possibility that a lunar satellite may have a stable orbit but that for some 
value g 0 of the argument of pericenter, dg/dtwill be zero so that dt/dgwill be 
undefined. In such an event the special values of 2g used to evaluate the coeffi- 
cients in the cosine series representation of 77 2 and dt/dg must be restricted to 
the region of admissible values, i.e., g 0 < g < (n - g Q ). 
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Table 3 

if, g; Elliptic Integrals and Harmonic Analysis (J 2 = 0) 


t 

(days) 

rj 2 

Elliptic 

Integrals 

V 2 

Harmonic 

Analysis 

■ 

g 

Elliptic 

Integrals 

(radi 

g 

Harmonic 

Analysis 

ans) 

■ 

0 

.84027778 

.84027778 

0 x 10' 5 

0 

0 

0 X 10“* 

20 

.83674774 

.83675070 

.296 

.22747356 

.2273639 

1.0966 

40 

.82691973 

.82692477 

.504 

.45445976 

.4543801 

.7966 

60 

.81298180 

.81298783 

.603 

.68513324 

.6850423 

.9094 

80 

.79820373 

.79821500 

1.127 

.92893740 

.9287134 

2.240 

100 

.78626473 

.78627455 

.982 

1.1954244 

1.1952296 

1.946 

120 

.78029580 

.78029692 

.112 

1.4844038 

1.4843190 

.848 

140 

.78191519 

.78191710 

.191 

1.7799258 

1.7799106 

.152 

160 

.79068037 

.79068695 

.658 

2.0605636 

2.0606739 

1.103 

180 

.80424132 

.80424327 

.195 

2.3167928 

2.3168499 

.571 

200 

.81911161 

.81911010 

.151 

2.5537699 

2.5537270 

.429 

220 

.83167986 

.83167979 

.007 

2.7820064 

2.7820046 

.018 

240 

.83906625 

.83906558 

.067 

3.0090287 

3.0090020 

.267 

260 

.83965270 

.83965478 

.208 

3.2366909 

3.2365252 

1.657 

280 

.83331311 

.83332190 

.879 

3.4638279 

3.4636120 

2.159 

300 

.82142913 

.82144059 

1.146 

3.6914969 

3.6913262 

1.707 

320 

.80670000 

.80671501 

1.501 

3.9263168 

3.9260582 

2.586 

340 

.79266806 

.79268607 

1.801 

4.1787443 

4.1783300 

4.143 

360 

.78290935 

.78291847 

.912 

4.4555575 

4.4552326 

3.249 

380 

.78002606 

.78002546 

.060 

4.7499516 

4.7497667 

1.909 

400 

.78480504 

.78480573 

.069 

5.0414831 

5.0413928 

.903 

420 

.79594658 

.79594596 

.062 

5.3120974 

5.3120934 

.040 

440 

.81050828 

.81050099 

.729 

5.5591054 

5.5589943 

1.111 

460 

.82483048 

.82482343 

.705 

5.7911936 

5.7910822 

1.314 

480 

.83552045 

.83551665 

.380 

6.0183377 

6.0182398 

.979 

500 

.84018034 

.84017939 

.895 

6.2456836 

6.2454964 

1.872 

520 

.83780104 

.83780875 

.771 

.19007482 

.1897622 

3.1262 

540 

.82889504 

.82891087 

1.583 

.41701462 

.4167136 

3.0102 
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* 


The harmonic analysis method makes it possible to include the oblateness 
of the moon when finding rj 2 and g as functions of time. Greater accuracy may 
be achieved by the retention of terms of the series after 8g or 8 \p. An advantage 
of this method is the ability to include the oblateness of the moon. This is espe- 
cially advantageous if the Orbiter is to be near the moon where the effect of the 
moon’s oblateness will be greater. 
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